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ABSTRACT 

Mid-infrared spectra of a few T Tauri stars (TTS) taken with the Infrared 
Spectrograph (IRS) on board the Spitzer Space Telescope show prominent nar- 
row emission features indicating silica (crystalline silicon dioxide). Silica is not a 
major constituent of the interstellar medium; therefore, any silica present in the 
circumstellar protoplanetary disks of TTS must be largely the result of process- 
ing of primitive dust material in the disks surrouding these stars. We model the 
silica emission features in our spectra using the opacities of various polymorphs of 
silica and their amorphous versions computed from earth-based laboratory mea- 
surements. This modeling indicates that the two polymorphs of silica, tridymite 
and cristobalite, which form at successively higher temperatures and low pres- 
sures, are the dominant forms of silica in the TTS of our sample. These high 
temperature, low pressure polymorphs of silica present in protoplanetary disks 
are consistent with a grain composed mostly of tridymite named Ada found in 
the cometary dust samples collected from the STARDUST mission to Comet 
81P/Wild 2. The silica in these protoplanetary disks may arise from incongru- 
ent melting of enstatite or from incongruent melting of amorphous pyroxene, the 
latter being analogous to the former. The high temperatures of ~ 1200K-1300K 
and rapid cooling required to crystallize tridymite or cristobalite set constraints 
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on the mechanisms that could have formed the sihca in these protoplanetary 
disks, suggestive of processing of these grains during the transient heating events 
hypothesized to create chondrules. 

Subject headings: circumstellar matter, infrared: stars, stars: pre-main-sequence, 
planetary systems: protoplanetary disks 



1. Introduction 



Silica (chemical name for the quartz group of minerals and a synonym for silicon diox- 
ide) is a major ingredient of the Earth's crust. The basic structural unit of crystalline silica 
is the tetrahedron, with the oxygen atoms at the corners and the silicon atom in the cen- 
ter. While in crystalline silica, structurally known as a tectosilicate, all oxygen atoms of 
a given tetrahedron are shared with the adjacent tetrahedra forming a three dimensional 
network structure, this does not hold for most other silicates. In silicates, Si04~ tetrahedra 
are coordinated by metal ions such as Mg^+ and Fe^"*"; examples of silicates are pyroxenes 
([Mg,Fe]Si03) and olivines ([Mg,Fe]2Si04). The spectra of both silica and sihcate dust have 
vibrational spectral features near 10 /im (due to the Si-0 stretching mode) and ~ 20 /xm (due 
to the 0-Si-O bending mode). For silica, the main features are at ~ 9, ~ 12.6, and ~ 20 /im, 
with a minor feature at ~ 16 /im for one of the polymorphs of silica, cristobalite, although 



Swainson fc Dovd (119931 ) note that the 16 /xm feature of /3-cristobalite is present at room 
temperature, but grows weaker with increasing temperature until it vanishes at ~ 525K. 
These are the features that we use to identify the presence of silica grains in our spectra of 
protoplanetary disks, which are disks of dust and gas thought to be raw material from which 
larger bodies, including planets, form. 

Silica has been used to dete rmine the temperature and pressure conditions at the time 
of forma tion in terre s trial rocks (Tuttle fc Bowerull958l ) and also in enstatite chondrite me- 



teorites (jsinns 1967 : 



Doddlll98lf). These conditions dictate how silica crystallizes into its 



various polymorphs (IHeaney et al.l Il994j : iHemley et al.l Il994l ) . Generally, for pressures of 1 
atm or lower, silica crystallizes as a-quartz at temperatures below 846K; as /3-quartz for tem- 
peratures between 846K and 1140K; as tridymite between 1140K and 1743K; as cristobalite 
between 1743K and 2000K; and as a liquid of silica composition above 2000K. For a phase 
diag ram illustrating thes e ranges of stability of silica polymorphs, see Figure 1 in Chapter 
1 of iHeaney et al.l (119941 ). At higher pressures and a wide range of temperatures, coesite, 
and then stishovite forms. All of the crystalline forms of silica are similar in that their bulk 
composition is Si02 and each of the oxygen atoms of every Si04~ tetrahedron is shared with 
the adjacent tetrahedron; however, the polymorphs differ in their arrangement of the Si04~ 
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tetrahedra in crystalline networks. If silica is not allowed sufficient time to crystallize into 
one of its polymorphs, it will form as amorphous Si02. In principle, one can estimate the 
formation conditions (temperature and pressure) of the silica dust in astronomical environ- 
ments based on the polymorph(s) of silica which can be inferred from their infrared (IR) 
spectra. 

Silica cannot be an abundant interstellar dust species sinc e the expected ~9 and 
~12.6/im bands are not seen in the ISM (see iLi fc Draind l2002l ). The presence of silica 
dust grains in the protoplanetary disk around the TTS Hen 3-600 A was indicated in its 
8-13 fim spectrum (IHonda et al.ll2003l ). 



In this work we report the detection of silica dust grains in the disks around five T Tauri 
stars from the prominent emission features characteristic of submicron silica grains. These 
emission features are generated in the optically thin uppermost layers of pr otoplanetary 



disks . We model the spectra obtained with the Infrared Spectr ograph (IRS; iHouck et al. 



20041 ) on board the Spitzer Space Telescope (I Werner et al.l 120041 ) with silica as a key dust 
component, with the aim to constrain the polymorph(s) of the silica dust in these disks and 
therefore the formation conditions (temperature and pressure history) of this dust. 



2. Data Reduction 
2.1. Observations 

The five T Tauri stars are described in Table 1. All were observed using both orders 
of the two low spectral resolution modules of the IRS, Short-Low (SL) and Long-Low (LL), 
using Staring Mode. For more information on the spectral modules and on Staring Mode, 



we refer the reader to lSargent et al.l (120061 ). These stars were chosen for our study of silica 
features because all five stars - IRXS J161410.6-230542, IRS 49, ROXs 42C, T51, and ZZ 
Tau - are SED Class II Young Stellar Objects; have prominent features at 9, 12.6, and 20 
fim wavelengths attributable to silica and often 16 /xm features partly attributable to silica; 
and have high quality spectra with high signal-to-noise ratios. 



2.2. Extraction and Calibration of Spectra 

We obtained basic calibrated data (BCD; fiat-fielded, stray-light-corrected, dark-current- 
subtracted) products from the Spitzer Science Center (SSC) for each of the targets in our 
sample from the S15.3.0 IRS data calibration pipeline for IRXS J161410. 6-230542 and from 
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the S14.0.0 IRS data calibration pip eline for IRS 4 9, RQ Xs 42C, T51, ZZ Tau. First, we 
identified bad pixels as described by IWatson et al.l (120071 ) in their explanatory supplement. 
In addition, we also identified as bad pixels a string of pixels with large values of e-/sec in 
the BCD data for the first nod of the first order SL observation of ZZ Tau from Campaign 
4. Some of the affected pixels were flagged as "radhit" detections (most likely a cosmic 
ray hit) in the associated bmask; we marked as bad pixels these and additional neighboring 
pixels that also appeared to be affected by the cosmic ray. As 12.5 /im is about the peak 
wavelength of a key diagnostic silica feature, fixing these pixels was critical for our purposes 
of analyzing dust emission features of silica. We correct fo r bad pixels by inter polation of 
neighboring pixels in the spectral direction, as described by IWatson et al.l (120071 ). First, we 
take the average of the pixel one pixel above and one pixel below. This is applied to all bad 
pixels in the bed data. If a given pixel is not fixed, we then take another iteration; this time, 
we take the average of the pixel one pixel above and two pixels below the bad pixel. This 
process is repeated, with subsequent iterations going two pixels above and one pixel below, 
then two pixels above and two below, then finally one more iteration going one pixel above 
and one below. 

The spectra presented in this paper were extracted fr om the BCD data u sing the Spec- 
tral Modelin g, Analysis, a n d Red uctio n Tool (SMART: iHigdon et al.l |200J) via methods 
described by Sargent et al. ( 2006 ) and Furlan et al. ( 20061 ). For all targets in our sample, 
we subtracted sky for SL and LL by subtracting the pixel-fixed BCD data for the same 
order but other nod position. We extracted a spectrum each DCE separately, then used the 
extracted spectrum from each DCE to compute, first, the average spectrum and, second, 
its uncertainties from the standard deviation of the mean. Wavelength calibration of our 
data, tapered column extraction for each order of SL and LL, and flux density calibration of 
spectra using Relati ve Spectral Response Functions (RSRFs) was achieved as described by 
Sargent et al.l (120061 ). The same high s pectral resolut i on tem plate of a Lac (Al V; M. Cohen 
2004, private communication) used by lSargent et al.l (120061 ) was used to generate RSRFs for 
SL and LL flux calibration for IRS 49, ROXs 42C, T51, and ZZ Tau, and for SL flux cali- 
bration for IRXS J161410.6-230542. For the LL flux calibration of IRXS J161410.6-230542, 
we applied LL second order, LL bonus order, and LL first order (from its short-wavelength 
end to 36 jxm wavelengt h) RSRFs generated from data of ^ Dra and the template for ^ Dra 
from lCohen et al.l J2003h : past 36 /xm, we applied a LL first order RSRF generated from data 
of Markarian 231 and the template for Markarian 231 (J. Marshall, private communication; 
Marshall et al1l2007l : lArmus et aDl2007[ l. 



We note the Humphreys-a and [Ne II] lines at 12.37 and 12.81 /xm, respectively, are 
located near the central wavelengths of the two narrow features consi tuting a double peak 
in the opacity profile of a-quartz (see discussion later in Section 3.1). ISargent et al.l (120061 ) 
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concluded a-quartz was present in the IRS spectrum of TW Hya based upon two narrow 
emission features located at the wavelengths of the two peaks of the a-quartz opacity pro- 
file used in t hat study. Upon de tailed inspection of the Short-High spectrum of TW Hya 
presented by lUchida et al.l (120041 ) , we find unresolved emission lines from Humphreys-a and 
[Ne H] and no underlying emission from the a-quartz double feature. We have looked at the 
Short-High spectra from Staring Mode observations of IRXS J161410.6-230542 (AOR ID# 
5453824) and ROXs 42C (AOR ID# 6369792) for these hues. The AORs had two DCEs per 
nod position and two nod positions; in these Short-High data, we found no emission from 
Humphreys-a or [Ne II]. 



2.3. Extinction Correction of spectra 



For all objects except the highly extinguished (^^=10. 7) IRS 49, we corre cted our 
RSRF-calibrated spectra in the same manner as described by lSargent et al.l (120061 ) . having 
adopted values of extinction at V-band, Ay, as listed in Table 1. For IRS 49, extinction 
correction was accomplished via the same method, and also using a lower rg.y assuming 
A-v / '^9.7=25, a value suggest ed for extinction in dense star forming regions like Ophiuchus 
(in which IRS 49 resides) bv lchiar et aP (12003). 



2.4. Mispointing 

We account for mispointing along the slits (the cross-dispersion direction) by adjusting 
the positions of our 3-5 pixel wide extraction boxes. The degree of mispointing is typically 
much less than 1 pixel (1.8 " in SL), and we adequately account for this by using the 
BCD flatfielding. Mispointing in the dispersion direction is harder to detect and causes a 
loss of signal at every wavelength. We assume that the flux density from one nod being 
systematically lower than that from the other nod was due to loss of signal from dispersion- 
direction mispointing. To correct for this signal loss, we multiplied the nod lower in flux 
density by a scalar to minimize the squared differences in flux at each wavelength in the 
order. For more about correction for mispointing, see the Appendix. 



2.5. Uncertainties 



Uncertainty in flux density was determined by the standard deviation from the mean 
for the scalar-multiplied spectra from all DCEs. For IRS 49, ROXs 42C, and T51, there 
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were only 2 DCEs (the 2 nod positions), so for these the standard deviation from the mean 
eq uals half the d i ffereii ce of the spectra at the 2 nod positions. Also, for the reasons listed 
by lSargent et al.l (120061 ) . we set relative uncertainties lower than 1% to 1%. 



2.6. Photosphere subtraction 



Three of the five objects whose spectra are analyzed in this study, IRXS J161410.6- 
230542, ROXs 42C, and ZZ Tau, have Spectral Energy Distributions (SEDs) indicative of 
very little ma ss in small dust grains in the innerniost d i sk regions. The S EDs of these 
three objects (ISilverstone et al.ll2006l : iMcClure et al.ll2008l : iFurlan et al.ll2006l ) have optical 
and JHK photometry that is consistent with stellar photospheric emission, with ZZ Tau 
photospheric even to L band (~3.5 /im). Also in each of these three cases, the excess 
above stellar photosphere as seen from the IRS spectrum at 5 /zm wavelength is very small. 
Therefore, photosphere subtraction is important for th ese spectra in order to isolate the 
emission from circumstellar dust (see ISargent et al.ll2006l ) at the short wavelength end of the 
IRS spectra. The low extinction to these three objects (see Table 1) gives us confidence that 
we can reasonably account for stellar photosphere emission in the IRS spectra. Therefore, 
for each we subtract a Planck function at the stellar effective temperature scaled to equal 
the flux at the dereddened J band flux of each object. At Spitzer IRS wavelengths (> 5 
/im), we are effectively subtracting the Rayleigh- Jeans tail of these Planck functions. The 
solid angles and temperatures used for the Planck functions representing stellar photosphere 
emission are listed in Table 1. For the stars T51 and IRS 49, there is significant excess at 
all IRS wavelengths, so no stellar continuum was subtracted. 



3. Analysis 



3.1. Silica Opacities 



Spitzer fc KleinmanI (jl96ll ) and IWenrich fc ChristensenI (119961 ) measured the complex 
dielectric functio n e(\) = ei{X) + ie2{\) for a-quartz at "room temperature" (~300K); 
Gervais fc Pirioul Jl975l ) obtained e{X) for a-quartz at 295 K, 505 K, 620 K, 785 K, and 825 K 
as well as for /?-quartz at 860 K and 975 K. 

To compare these dielectric functions, we calculate, for the sake of ill ustration, the opac 
ities for g-quart z with the sanie gra in shape distributions of either CDE (IBohren fc Huffman 



19831) or CDE2 (IFabian et al.ll200lf). using the dielectric functi ons from lSpitzer &: Kleinman 
( jl961 ) for quartz at ~300K, Wenrich fc Christensen (1996) for quartz at ~300K, and 
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Gervais fc Pirioul (Il975l ) for quartz at 295 K and taking the "| — ^-approximation" for the 
ordinary and extraordinary rays (see § A in ISargent et al.l 120061 ). We find that the opacities 
computed from the three sets of dielectric functions are nearly identical. This implies that 
we can have confidence not only in the optical properties of a-quartz at r oom temperature, 
but al so in the dielectric functions of a-quartz and /3-quartz reported by iGervais fc Piriou 
(I1975I ) at other temperatures. 



We show opacity curves for CPE and CDE2 shape distributions using the dielectric 
function from lGervais fc Pirioul (1l975l ) for g-quartz at 295 K in Figure 1. We take the mass 
density of a-quartz to be 2.65 gcm~^ (IHeaney et al.l 1 19941 ). In the literature, the CDE and 
CDE2 grain shape distributions have been widely considered to explore the effects of dust 
shapes on the opacity profiles. Both CDE and CDE2 assume that the size of the grains under 
consideration are in the Rayleigh limit; i.e., 27ra/A << 1. Thus, the grain sizes for CDE 
and CDE2 shape distribution s used in our models of m id-infrared emission are submicron. 
The CDE shape distribution (IBohren &: Huffmanlll983l ) assumes a distribution of ellipsoidal 
shapes in which each shape, defined by the ax ial ratios of the el lipsoids, is statistically 
weighted equally. The CDE2 shape distribution (iFabian et al.ll200ll ) assumes a distribution 
of ellipsoidal shapes in which ellipsoids that are almost spherical are weighted more heavily 
than more extreme (highly elongated along one or two axes or highly flattened along the 
other axes or both) shapes; the most extreme shapes in CDE2 are given zero weight. CDE is 
therefore more heavily weighted toward very nonspherical ellipsoids than CDE2. Our results 
confirm what others have s hown in the literat ure (e.g. see the opacities of spheres, CDE, 
and CDE2 for forsterite by iFabian et al.ll200ll ) - shape distributions of ellipsoids that more 
heavily weight toward very nonspherical ellipsoids with respect to spheres result in opacity 
profiles with features centered at increasingly longer wavelengths. 

To approximate the effects of quartz of a range of temperatures, in Figure 2 we show the 
/5-quartz opacity profile obtained by averaging the opacity curv es of /3-quartz of 97 5 K and 
860 K. We take the mass density of /3-quartz to be 2.53 gcm~^ (IHeaney et al.lll994l ). Simi- 
larly, we obtain the high-temperature a-quartz profile from averaging the opacity of a-quartz 
at 825 K and 785 K, and the lower-temperature a-quartz profile from averaging the opacity 
of a-quartz of 620 K and 505 K, all with a CDE shape distribution. Note the progression 
of the 12.6 /im feature complex from double-peaked to single-peaked as the quartz temper- 
ature increases (see Fig. 2). We assume that a given polymorph of silica, once formed, can 
exist stably at temperatures and pressures within its stability field on the phase diagram or 
metastably outside the stability field at lower temperatures if the physical conditions change 
relatively rapidly. However, we assume a polymorph of silica becomes another polymorph 
if subjected to temperatures higher than allowed by its stability field in the silica phase 
diagram. 
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Infra red transmission spectr oscopy by iPlyusnina et al.l (119701 ) and absorption measure- 
ments by IWilliams et al.l (Il993l ) for coesite and stishovite rule out these higher pressure 
polymorphs o f silica in our sample of TTS disks. For the i n terme diate pressure polymorph, 
coesite, both IPlyusnina et al.l (Il970l ) and IWilliams et al. Jl993h show four features with 
widths similar to that of the 12.6 /im feature between 14 and 23 /im, but all with strengths 
greater than that of the 12.6 /xm feature - one feature near 14 /xm, a pair of slightly stronger 
features at ~ 16-17 /im, and the strongest of the four features at ~ 20-23 /zm. The spectra 
with prominent silica features do often have a feature at ~ 16-17 /im, but it is single, and, 
in addition, such spectra do not have a narrow feature at 14 /im of greater prominence than 
the 12.6 Atm feat u re. T his ru les out coesite. Stishov ite can be more easily ruled out, as both 
Plyusnina et al.l (Il970l ) and IWilliams et al.l (Il993l ) show stishovite to lack a 9 fim feature, 
while the spectra of the five silica exemplars analyzed in this study have prominent 9 fim 
features. 

To represent the two higher temperature, low pressure polymorphs of silica, tridymite 
and cristobalite, we use the opacity obtained from the transmission me asurements of subm i- 
cron annealed silica grains embedded in Potassium Bromide, KBr, by iFabian et al.l (120001 ). 
This silica was formed by heating amorphous Si02 grains at 1220 K for 5 hours, and was 
found to be mostly cristobalit e but also partial l y tridy mite . This is supp o rted b y the fact that 
the cristobahte presented by IPlyusnina et al.l (Il970l ) and iRokita et al.l (119981 ) have a 16 /xm 
feature of prominence similar to but slightly less than that of the 12.6 /xm feature, whereas 
for tridymite the 16 A tm fea ture is almost nonexistent (see the t ridymite abso r bance shown 
by iHofmeister et al.l (Il992l ) ) . The annealed silica presented by iFabian et al. J2000h has a 



16jxm feature of prominence closer to that of "low cristobalite" as me asured bvlRokita et al. 
(1998). We also note that the opacity profile of the silica annealed by Fabian et al. ( 200ol ) at 



1220 K for 5 hours bears a close resemblance to the IR spectrum of the pr oducts of anneal 



i ng fo r 120 hours at 1300 K of silica formed from iron-bearing silicates by iHallenbeck et al. 



(119981 ). We present in Figure 3 the opacity of the annealed silica used in our models. 



For amorphous S i 02, the dielectric fu i iction s at r aid-IR wavel e ngths were measured by 
Henning fc Mutschkd (119971 ) , ISteyer et al.l (119741 ) , and iKoike et al.l (Il989l). In our dust mod- 



eling w e adopt the dielectric functions of amorphous Si02 at 300 K oflHenning fc Mutschke 
(119971 ). In addition, we consider obsidian, as measured by iKoike et al.l (119891 ). To compute 
the opacities of obsidian a nd amorphous Si Og, we take their ma ss densities to be 2.384 g cm~'^ 



( Koike &: Hasegawalll987l ) and 2.21 g cm" -3 JPabian et allboooh . respectively. We plot in Fig- 
ure 3 also the opacities of obsidian and amorphous Si02 in the CDE shape distribution. The 
shape distribution affects the opacities of obsidian and amorphous Si02 less than it does for 
a- and /5-quartz; for this reason, we only use the opacities of obsidian and amorphous Si02 
with a CDE shape distribution, and we do not plot the opacities of these grains in the CDE2 
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shape distributions. 



3.2. The KBr Effect 



Embedding smaU grains in KBr is commonly used to measure the grain opacity in the 
IR. This technique has hmitations - optical constants are not directly measured, so it is 
not possible to explore grain shapes, porosity, or size. The embedded grains may clump 
together, reducing the contrast in sharp features. Moreover, the index of refraction of KBr 
(~ 1.5) can distort strong resonance features in the spectrum from how they would appear 
in a spectrum of the same grains in a vacuum (i.e., in a protoplanetary disk). 

Annealed silica and crystalline pyroxene are the only grain components in our models 
which depend on a KBr-embedded measurement. We model the KBr effect on /3-quartz 
(see Figure 4) to illustrate the probable effects of KBr on our annealed silica opacities. We 
compute the opacities of solid, small /3-quartz grains in a KBr medium and in vacuum for 
both CDE and C DE2 shape distributions . We use the optical constants of /3-quartz at 
975K computed by lGervais &: Pirioul (Il975l ) as the optical properties of silica grains that ar e 
embedded in a KBr medium whose 6-28 nm optical constants we take from iHeatonI (jl97ll ). 
As shown in Figure 4, for the CDE shape distribution, the peak wavelengths of the 9, 12.6, 
and 20 /im features are little altered by embedding in KBr. The peak wavelength of the 9 fim 
feature for /3-quartz in the CDE2 shape distribution embedded in KBr is at a slightly longer 
wavelength than that for /3-quartz grains in the CDE2 distribution in vacuum. This also 
holds for the 20 /im feature. The peak wavelength of the 12.6 fim and 16 fim features are 
much less affected by either being in KBr or with a different shape distribution. 

Note in Figure 4 that if the shape distribution for the dust in KBr is already at an 
extreme (e.g. CDE), models with the same shape distribution but in vacuum will have 
spectral features not much different (in width and central wavelength) from those in KBr. 
We note that all opaci t ies of /9-quartz are somewhat higher for dust in a KBr medium, also 
found by IJaeger et al.l (119981 ) . 



3.3. Other Dust Opacities 



For submicron amo rphous silicates, we compute opacities for the CDE2 shape distribu- 



tion (Fabian et al.ll200ll ). We take the dielectric functions of MgFeSi04 of iDorschner et al. 
(119951 ) for amorpho us olivine. For amorph ous pyroxene we use the dielectric functions of 
Mgo.rFeo.sSiOa from IDorschner et al.l (119951 ). We take the mass densities of both olivine and 
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pyroxene to be 3.3 gem 



To account for the effects of grain growth, we consider poro us dust. Porous material i s 
simulated by applying the Bruggeman effective medium theory (IBohren fc Huffman! Il983l ). 
assuming the porous grains are 60% vacuum by volume. The effective di electric functions for 
porou s amorphous pyroxene and olivine are then computed. Mie theory (IBohren fc Huffman 
19831 ) is t hen used to compute the opacities of spherical porous grains of radii of 5 /im. As 
noted by ISargent et al.l (120061 ) . the opacity of 5 /xm-radius, 60%-porous ar aorphous silicate 



grains is similar to that for 2 /im-radius solid dust, the size of dust used by iBouwman et al. 
(I2OOII ) to investigate grain growth in Herbig Ae/Be disks. 



We also include pyroxene (crystalline) and forsterite (of whic h the density is also 3 .3 g cm~ 
in our analysis. For pyroxene, we use the opacity of "En90" (jChihara et al.ll2002l ). With 
a stoichiometric composition of Mgo.gFeo.iSiOs, "En90" is quite effective in fitting the py- 
roxene features at 9.3, 10.5, 11.2, and 11.6 /xm within the lO/ii n complex of FN Ta u (see 
Sargent et al.ll2006l ). For forsterite, we use the optical constants of lSogawa et al.l (120061 ) which 
are similar to that of IServoin fc Pirioul (Il973r ). to compute the opacity for a distribution of 
ellipsoidal shapes. 



Using the formalism of iBohren &: Huffman! (!l983! ). we have constructed an opacity pro- 
file for a distribution of forsterite grains in which all ellipsoidal axial ratios are given an 
equal weight, as in CDE. However, we exclude the most extreme ellipsoidal axial ratios in 
the following manner. 

Our shape distribution, which we call "tCDE" (truncated Continuous Distribution 
of Ellipsoids), is that we confine Li and L2 (which are the Lj parameters described by 
Bohren fc Huffman! (119831 )) within the triangle in the Li-L2-space specified by vertices {Li = 
0.1, L2 = 0.05), (0.1, 0.895), and (0.99, 0.005). For CDE, all shapes are allowed: the triangle 
specified by the verticies (0,0), (0,1), and (1,0). The opacity profile of forsterite for CDE 
has a feature peaking at 11.3 /xm which is close to the feature seen in astronomical data (but 
often peaking at a slightly longer wavelength). The opacity profile for tCDE "rounds" this 
peak somewhat, pushing the peak toward slightly shorter wavelengths, and this results in 
marginally better fits to astronomical spectra. 



We will discuss "tCDE" at a greater length in !Sargent et al.! (12008! ) . but here we justify 
this choice by showing in Figure 5 our model fits to ROXs 42C with forsterite grains having 
either a CDE shape distribution ("CDE" model) or a tCDE distribution ("tCDE" model). 
The former (CDE) has psr degree of freedom (d.o.f.) of ~ 5.1, while for the latter (tCDE) 
X^/d.o.f. ~ 4.3 which is somewhat better; the tCDE model does a slightly better job in 
fitting the 11 and 23 /xm complexes of the spectrum of ROXs 42C but is overall very similar 
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to CDE. ISargent et al.l (120061 ) found great similarity between the opacity of forsterite and 
silica grains with a CDE shape distribution and that of porous dust with a porosity of 60%. 
Thus the effects of shape distribution and porosity cannot be clearly distinguished in the 
spectra. 

No hint of Polycyclic Aromatic Hydrocarbon (PAH) emission is seen in our spectra, 
except possibly ZZ Tau at 6.2 fim. Especially telling is the lack of the ubiquitous 6.2, 
7.7/im PAH features, which are clearly distinguishable from silicate or silica emission; in ZZ 
Tau, no 11.3 /im PAH emission is evident. 

We note that the opa city profiles of so r ae sol ids are similar to that of silica, but nev- 
ertheless can be ruled out. iKimura fc NuthI (120071 ) measured the transmittance of H2Si204 
and Si203. They found in some of their spectra the features peak at wavelengths very 
close to those of silica (~9, 12.6, and 21.5 /zm); however, many of their spectra also show 
a narrow feature at 11.36 /zm with a width similar to their 12.6 /xm feature but of greater 
prominence. This feature is not seen in the spectrum of ZZ Tau. In the spectra of IRXS 
J161410. 6-230542 and ROXs 42C there is a feature near 11.1 /im, but as will be shown in our 
model fits to their spectra, this peak can be attributed to forsterite. Also, the opacity curves 
cor nputed using the optic al constants of Na-bearing amorphous aluminosilicates determined 
by iMutschke et al.l (119981 ) exhibit features near 9 fim and 21 fim similar to those of silica of 
various versions; however, the Na-bearing aluminosilicates lack the distinctive 12.6 /xm fea- 
ture which is seen in the emission spectra of the 5 T Tauri stars and in the opacity spectra 
of all versions of silica. 



3.4. Models 



We model the spectra of each of five T Tauri stars as a sum of featureless continuum 
emission (from optically thick disk midplane, iron or amorphous carbon grains, and very large 
>10 nm silicate grains) and optically thin emission from dust grains with strong infrared 
resonances in the disk atmosphere (see lSargent et al.ll2006l ). It is not our goal to model self- 
consistently the entire spectral energy distribution from the nea r -IR to subr nm; such would 
requir e radiative transfer models such as those of ICalvet et al.l (119921 ) and iD'Alessio et al. 
J200lh . 



We reject silicate self-absorption, a mix of emission and absorption by similar kinds 
of dust but with cooler dust in front of warmer dust, as an explanation of the 9 and 12.6 
fim features we see in our spectra. Self-absorption in the 10 fim silicate complex would 
give apparent emission features in the 8 and 12 /xm wings of the 10 /xm complex; however. 
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such would require the optical depth of the cool amorphous silicate dust to exceed 1. Since 
amorphous silicate dust in the ISM gives Ay jr^ -j of 18 or 25 (see §2.3), this would mean the 
intervening dust has Ay of 18 to 25. Such dust would not be at the disk surface; otherwise 
it would heat up to the same temperature as surrounding dust and not give rise to an 
absorption feature. Such dust would have to be further away from the disk and star, and 
therefore would also be in front of the star in the line-of-sight to the star, and we would 
therefore find Ay of 18 to 25 to all stars with what appear to be prominent 9 and 12.6 /im 
silica features. All of our objects have Ay less than 18, the minimum Ay required to result 
in self-absorption. 

We approximate the IR emission from dust wit h a large range of te mperatures by two 
temperatures. A similar approach was employed by iKastner et al.l (120061 ) to model the IRS 
emissio n spectra of the d ust around the hypergiant R66 in the Large Magellanic Cloud 
and by I Chen et al.l (120061 ) to model that of five debris disks (see especially the two dust 
temperatures used in the model for 77 Tel). We assume that, for a given protoplanetary 
disk, all dust grains, regardless of composition, within the same arbitrarily small volume of 
the part of the disk atmos phere giving rise to opticall y thin silicate and silica emission are 
at the same temperature. iKamp fc DuUemondl (120041 ) find in their disk models that, in a 
layer above the layer giving rise to dust emission features, the gas and dust temperatures are 
within 10% of each other. As gas density increases closer towards the disk midplane, the gas 
and dust temperatures in the layer giving rise to dust emission features should be even closer 
to each other. The dense gas effectively should equalize the temperatures of dust grains of 
different optical properties (resulting from, e.g., different compositions, different sizes). Our 
two-temperature models therefore model dust composition for two disk regions: the warmer 
inner disk regions and the cooler outer disk regions. 

Our model is a sum of blackbodies at two temperatures plus optically thin emission 
from dust at those temperatures. Mathematically, these models are given by 



mod 



(A) + 



(1) 



where Fy{X)^°'^ is the model flux density; By{\^ T) is the Planck function at temperature T 
and wavelength A; Tc is the temperature of "cool" dust and is the temperature of "warm" 
dust; fic (^wj) is the solid angle of a blackbody representing the continuum emission from 



cool (warm) dust; a^j 
where rrici 



), the mass weight for cool (warm) dust, equals mc,i/d^ 



m^j) is the mass of dust species i at Tc (T^) and d is the distance to the T 
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Tauri star (see Table 1 for the stellar distances). We specify a grid of temperature pairs, 
Tc and T^, not allowing either of the temperatures to exceed 1401 K, as this temperature 
is intermedia te between the sta bility limits of forsterite and enstatite of 1430 K and 1370 K, 
respectively (jPosch et al.l 120071 ). For a given temperature pair, our model equation is linear 
in solid angles and optically thin dust mass weights. We minimize 



with respect to each of the solid angles and mass weights, where Fy(Afe)"'^ is the observed 
flux density at the k-th waveband and AFi^lX^y^^ is the uncertainty in the observed flux 
density at the k-th waveband. If any of the mass weights computed from an iteration of 
minimization are negative, the component with the most negative integrated flux from 7.7-23 
/xm is set to zero, and we then minimize again. This process is iterated until there are no 
components with negative mass weights. The solution found in the iteration for which no 
components are negative is the best fit for a given temperature pair. The process is carried 
out for another temperature pair. The global best fit for a spectrum is at the temperature 
pair that minimizes P^r degree of freedom. In all of our best fits, the silica component 
at both low and high temperatures was positive. The number of degrees of freedom is equal 
to the number of data points minus four (since we consider two dust temperatures and two 
solid angles) minus the number of mass weights at both temperatures altogether. 

Figure 5 is a sample fit of a model with annealed silica to the 7.7-37/im IRS spectrum 
of ROXs 42C (all other models, including the other models of ROXs 42C, discussed or shown 
in this paper are fit over 7.7-23 yum). It is seen that our two-temperature model fits the 
spectrum of ROXs 42C over almost the entire wavelength span of the IRS. What interests 
us most in this study are the spectral features diagnostic of silica, the 9, 12.6, 16, and 20 /im 
features. For this reason, we minimize over the smallest range of wavelengths sufficient to 
adequately cover these features, 7.7 to 23 /im. We determine a very conservative estimate of 
the uncertainty for each nonzero dust component and solid angle in the model as the amount 
that resulted from increasing per degree of freedom by ~ 1 from the best-fit value. The 
models shown in all figures after Figure 5 and whose parameters are summarized in Table 2 
are constructed by minimizing over 7. 7-23 fim. 

We obtained models of quartz in each of the CDE and CDE2 shape distributions (where 
the same shape distribution is used for quartz at the low temperature, Tc, as for quartz at 
the high temperature, T^). We used the opacity of the 295 K a-quartz for the "cool" silica 
component. For the opacity of the "warm" silica component, we use the /5-quartz, high- 
temperature a-quartz, and lower-temperature a-quartz opacity curves described in §3.1 in 



-14- 



reference to Figure 2. The opacities computed from the dielectric constants of a-quartz 
at 505 K and 620 K were averaged to give an "average 562.5 K a-quartz opacity"; 785 K 
and 825 K a-quartz opacities were averaged to give an "average 805 K a-quartz opacity" ; 
and 860 K and 975 K /5-quartz opacities were averaged to give an "average 917.5 K /5-quartz 
opacity". The model with the lowest difference between T^j and "average a- or /?-quartz 
opacity temperature" (562.5 K, 805 K, or 917.5 K) is the best fit model for quartz grains in 
the shape distribution under consideration (only CDE and CDE2 were considered) . A model 
is rejected if is greater than the highest temperature allowed according to the stability 
field for a given silica polymorph (for silica polymorph stability fields, see Section 1). 



Results 



4.1. ROXs 42C 



ROXs 42C was identified by iBouvier fc Appenzellerl (119921 ) as the easternmost stellar 
source within the positi o nal error of the X-ray source ROX 42 found by the Einstein satellite 
( jMontmerle et al.lll983l ). iMathieu et al.l (119891) f o und it to have a double-lined spectroscopic 



binary with a 36 day period, while 
separated by ~ 0.16": 



GhezetaL 



(119931 ) found it also to be a visual binary 



Jensen fc Mathieul (119971 ) conclude that ROXs 42 C is a hierarchical 



triple - a spectroscopic binary orbited by single star further out. iRatzka et al.l (120051 ) agree 
with the assessment that ROXs 42C is a hierarchical triple, but report that the visual binary 
is now separated by 0.277". Assuming a distance to Ophiuchus of 140pc (see Table 1), this 
implies a se paration of ^ 39 A U between the double-lined spectroscopic binary and the 
single star. IJensen et al.l (119961 ) calculate the separation of the double-lined spectroscopic 
binary to be 0.27 AU. We note the similarity between ROXs 42C and Hen 3-600, itself a 
spectroscopic binary wi th another star at 70 AU orbiting the binary. We also note that 



(Honda et al. 


20 


03; 


Uchida et al. 


2004) 


is unknown, as 


Bouvier & Appenzeller 



neither in absorption or emission, but such that the photospheric line was filled in by 1.6 A 
of emission. However, because of the multiple stars in this system, and because it is not 
known where the disk(s) is (are) located, we do not know whether active accretion could be 
occurring currently onto one star and its spectral signature is obscured, or whether there is 
no active accretion in this system. 

The IRS spectrum of ROXs 42C has prominent narrow features at 9 and 12.6 fim from 
silica, in addition to features around 11, 23, 28, and 33.5 fim identified with forsterite and 
pyroxene. We show the model fits to ROXs 42C using each of annealed silica, obsidian. 
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amorphous Si02, and quartz grains in the CDE and CDE2 shape distributions in Figure 6. 

We rule out the models using the amorphous forms of silica - obsidian and amorphous 
Si02 - because the 12.3-12.4 //m features of these two sihcas are of much less strength with 
respect to their 9 //m features than for the crystalline forms of sihca. This translates into 
being able to fit the 9 and 20 fj,m features somewhat reasonably but not having enough 
strength in the 12.3-12.4 //m feature to fit the data. In addition, the 12.3-12.4 feature 
from the amorphous versions of silica is centered shortward of the feature in the data by ~ 
0.2 fim and is noticeably wider, with a Full Width Half Maximum (FWHM) of ~ 1 /im, than 
that in the data (with a FWHM ~ 0.5 fim.). The pcr degree of freedom for the model 
using obsidian is 5.4; for the model using amorphous Si02, it is 5.2. 

Next, consider a-quartz, /3-quartz, tridymite and cristobalite (tridymite and cristobalite 
being represented together by annealed silica). For quartz in the CDE2 shape distribution, 
the lowest difference between model and "average opacity temperature" (temperature at 
which quartz optical constants were measured) happened for /3-quartz. A per degree of 
freedom of 4.5 was obtained for this model. The fit to the 9 //m feature is good, coming 
from a combination of /3-quartz and pyroxene. Also, the model is within the error bars of 
a httle over half the data points between 20.5 and 22 /xm. However, the spectrum of ROXs 
42C lacks the double feature at 12.6 /im in the model, and the model does not have a strong 
16 fim feature like ROXs 42C. Forsterite gives a 16 //m feature, but forsterite is insufficient 
by itself to fit the 16 jj,m feature in ROXs 42C. 

A per degree of freedom of 5.0 was obtained for the model using quartz in the CDE 
shape distribution. The fit to the 12.6 /im feature is improved due to the feature being 
single-peaked in higher-temperature a-quartz. However, this model provides a poor fit to 
20.5 - 22.0 fim. The model feature at 9 fim peaks slightly longward of the feature in the data. 
A shape distribution intermediate between CDE and CDE2 could fit the 9 fim feature, but 
would still provide a 20 /xm feature that would overshoot the ROXs 42C spectrum between 
20.5 and 22.0 iim. Furthermore, the model does not have a 16 iim feature sharp enough to fit 
the 16 /im feature in the data. We believe this definitely rules out /3-quartz as a candidate 
to explain the sihca features of ROXs 42C; a-quartz at temperatures of 505K-620K is a 
possibility to explain the features, but we believe it to be of lower probability. 

This leaves annealed sihca. We show the model using annealed silica at low and high 
model temperatures in Figure 6. A per degree of freedom of 3.2 was obtained for this 
model, the lowest for any of the models of ROXs 42C. The fit to the 9 //m, 12.6 //m, and 

16 nm features and the 20.5-22.0 fim. region is the best of all models fit to ROXs 42C, and 
this is refiected in the x"^ pcr degree of freedom (Table 2 and Figure 7). We note the opacity 
of annealed silica that we use comes from KBr pellet transmission measurements, but, as 
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discussed in Section 3.2, we expect that the opacity curve for the annealed sihca grains in 
KBr would be very similar in terms of spectral feature central wavelength and width to the 
opacity in vacuum. We conclude that annealed silica best matches the spectrum of ROXs 
42C. 



4.2. ZZ Tau 



Simon et al.l (Il996l ) found ZZ Tau to be a binary system with a separation of stars 
that increased from 0.0338" in 1994 to 0.0420" in 1996. S ubsequent measurements of the 
positions of the stars in the biriary let ISchaefer et al.l (120061 ) solve for the orbital parameters 
of this system. ISchaefer et al.l (120061 ) com pute the semi-rn a ior a xis expressed as an angle 
of 0.061", implying a separation of 8.5 AU. iPrato fc Simon! (119971 ) report EWHa from this 
system of 15 A and suggest that this system is either a hybrid between a weak-lined and a 
classical T Tauri Star (the disk around one star actively accreting, the disk around the other 
not) or that one of the two stars in this system has an inner gap (a "transitional disk"). 

ZZ Tau was also fit by amorphous Si02, obsidian, quartz in both CDE and CDE2 shape 
distributions, and annealed silica. For this object and for the three other objects, we only 
show the best fit. As with ROXs 42C, amorphous Si02 and obsidian gave features much 
too weak, too wide, and centered at the wrong wavelength to fit the 12.6 fim feature of 
ZZ Tau. The fits using quartz in both shape distributions were rejected because the high 
temperature in the model was beyond the range of stability for /?-quartz. The fit to ZZ Tau 
using annealed silica (Figure 8, top) fit quite well at 9, 12.6, 16, and 20 /im, and it also 
proved to be the best fit in terms of P^r degree of freedom, which amounts to 3.8 (see 
Table 2 and Figure 7). 



4.3. IRXS J161410.6-230542 



Of the 7 sources reported bv lMetchevI (120061 ) to be located within 10"of IRXS J161410.6- 
230542, only one was confirmed to be associated with this star. This companion is .222" 
separated from the primary and is 0.21 magnitudes fainter at K-band (jMetchevll2006l ). At 
an assumed distance to t his object of 145pc (T able 1), this implies a separation of 32 AU 
between the components. IPreibisch et aD Jl998h report EWHa from IRXS J161410.6-230542 
of 0.96 A, whi c h, wh en combined with the KO spectral type assigned to this system by 
Pascucci et al.l (120071 ). suggests little or no accretion is occurring from the disk to the star 
in the IRXS J161410.6-230542 system. 
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Annealed silica provided the best fit to the Spitzer IRS spectrum of IRXS J161410.6- 
230542. As seen at the bottom of Figure 8, the fits to the 9, 12.6, 16, and 20 /im features are 
all adequate with our model. The P^^' degree of freedom is 6.6 for this model. Of all the 
versions of silica discussed in this paper, only annealed silica has a single, isolated peak at 
16 /im. The 16 /im feature in our model is coming equally from forsterite and cold annealed 
silica, and the model fits the narrow (~ 1 /im FWHM) peaks in the spectrum fairly well at 
11 and 19 /im, which are features associated with forsterite, meaning the fit to the 16 /xm 
feature would be worse without the contribution from annealed silica. 

We note that the curious trapezoidal shape of the 12.6 /im feature does not perfectly 
match the rounded shape of the 12.6 /im annealed silica feature in our model, but we also 
note the size of the errorbars in the data at the top of the 12.6 /im feature do not preclude 
a more rounded shape like the feature belonging to annealed silica. The 12.6 /im complex of 
a-quartz does evolve from double-peaked at 295K to somewhat trapezoidal and eventually 
single-peaked at progressively higher temperatures of 505K, 620K, 785K, and 825K; the 
feature is single-peaked and centered at 12.7 /im for /3-quartz. As opposed to annealed 
silica, a- and /3-quartz in the CDE shape distribution do not fit the 20 /im feature, and 
both quartzes in the CDE2 shape distribution do not fit the 9 /xm feature (see discussions in 
sections 4.1 and 4.2 on ROXs 42C and ZZ Tau, respectively). Further, models with annealed 
sihca fit the 16 /im features of IRXS J161410.6-230542 and ROXs 42C noticeably better 
than does quartz in either the CDE or CDE2 shape distributions. 



4.4. IRS 49 



IRS 49, also kno wn as WLY 2-49, has been searched for m ultiplicity but has n ot been 
found to be multiple (IBarsony et al.ll2003l : iRatzka et al.ll2005l ). iRatzka et al.l (120051 ) estab- 
lished upper limits on the relative brightness at K band of any unseen companion around 
IRS 49 of 0.04 for 0.15" separation and 0.02 for 0.50" separation. ICatti et all J2006h report 
significant accretion of M = 1.3 x 10^^ MqUt'^ from Paschen-/3 and Brackett-7 emission. 

IRS 49 is in the Ophiuchus star-forming region; it has a sizeable extinction of Ay of 
10.7 (Table 1). How extinction at visible wavel engths transla tes to extinction at mid-infrared 
wavelengths is still uncertain (see discussion by lDraind (120031 ) on values of Ay/rg 7 determined 
from various studies); for this reason, we present fits to the spectrum of IRS 49 dereddened 
assuming Ay/T9.7=18 (Figure 9, top, brown points and solid black line; per degree of 
freedom of 5.5) and assuming Ay /rg. 7=25 (Figure 9, top, orange points and solid violet line; 



per degree of freedom of 6.3). As discussed by ISargent et al.l (120061 ) with regard to the 



dust model of V410 Anon 13, greater extinction correction of at mid-infrared wavelengths 
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results in spectral profiles that look more like the interstellar medium profile; namely, more 
like amorphous silicates. The effect of extinction correction on the silica features is most 
important for the 9 and 20 /xm silica features, but, as can be seen from Table 2 and Figure 7, 
per degree of freedom is lowest for both extinction corrections using models with annealed 
silica. 



4.5. T51 



Lastly, we consider T51 . Also known as Sz 41, this object was found to be double by 
Reipurth &: Zinneckerl (119931 ). Another star, named Sz 41C, had a lso been thou ght to be 
associate d with Sz 41, but w as shown to be a background star by IWalterl (119921 ). This is 
noted b ylCorreia et al.l (120061) . whose observations show the Sz 41 binary to be separated by 
1.974". iGuenther et al. ( 2007) ass ign EWHa of -2.0A, i.e., absorption. Spectroastrometry 



conducted by iTakami et al.l (120031 ) determine that the Ha profiles from the two components 
of the Sz 41 binary are roughly equivalent, suggesting there is negligible accretion in the 
disk(s) of the Sz 41 system. 

We fit models with different silicas to T51. The spectrum of T51 includes a very promi- 
nent peak at 9 fim, a triangular 20 fim feature, a modest 12.6 fim feature, and a slight 16 
/im feature. Once again, the model with annealed silica fit the dust excess spectrum the 
best (see Figure 9, bottom; P^i' degree of freedom of 3.6). However, we note that of all 
objects whose spectra were analyzed in this study, T51 provides the weakest tests discrimi- 
nating between the various versions of silica (see Figure 7). The models using obsidian and 
amorphous Si02 are not much worse than the model using annealed silica. 



5. Discussion 

5.1. Cristobalite Dominance 

All the models (see Table 2) indicated the presence of both amorphous and crystalline 
silicates and grains grown to larger sizes as well as silica, suggesting production of silica 
accompanies both dust processing and grain growth. Models using the opacity of annealed 
silica provided the best fit of all models for all five objects whose spectra we analyzed. From 
our discussion of the KBr effect in Section 3.2, we believe the opacity of the same silica 
grains in vacuum instead of KBr would be very similar; any difference could be compen- 
sated by assuming more extreme ellipsoidal shapes (as in CDE) . We conclude that the silica 
polymorph giving rise to the silica features we see in the spectra of ROXs 42C, ZZ Tau, 
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IRXS J161410. 6-230542, IRS 49, and T51 is predominantly cristobalite becaus e cristobalite 



was th e dominant polymorph of silica present in the annealed silica sample of iFabian et al. 
( 120001 ). although an admixture of tridymite is possible. 



5.2. Silica in the Solar System 



Dodd 


( 


1981 


) rei 


also 


Binns Il967[ 



more specifically, they were found in ensta tite chondritic me te orites, 



which are characterized by calcium-aluminum inclusions (CAIs; iGuan et al.ll2000l ). iDodd 
( ll98ll ) also summarizes findings on eucrite meteorites. Eucrites are achondrites, which means 
they lack chondrules. Chondrules are spheroidal inclusions in chondritic meteorites, are 
submil limeter to cen timeter in diameter, and have the appearance of having formed as molten 
drops (jNortorul2002l ). The findings show the volume percentage of silica (as quartz, tridymite, 
and cristobalite polymo rphs) in these eucrites to be between and 4 %. In the review by 
Brearley fc Joned (119981 ). it is mentioned that less than 2% of chondrules are silica-bearing 
chondrules, chondrules that are up to 40% by volume silica, and they have been observed in 
a number of ordinary chondrites. 



Silica has not 



jeen widely reported as detected in the infrared spe ctra of come t co- 



mae (IHannerl I2OO3I : iHanner fc Bradley! |2004| ). Neither Comet Kohoutek (lMerrilllll974l ) nor 
Comet Hale-Bopp (IHarker et al.ll2002l ) showed any hint of the 9 and 12.6 /im silica features. 
Spitzer IRS spectra of the comet nucleus ejecta from the Deep Impact mission to Comet 
Tempel 1 also show no hint of sili ca features, and models of emission from this dusty eject 
do not require silica for a good fit (ILisse et al.l 120071 ). iBouwman et al.l (1200 ll ) reported silica 
present based on the spectrum of Comet Halley; however, the data is not of high quality. 



being coarse and noisy, so the detection is tentative at best. As discussed by iHarker et al. 



(I2OO2I ). magnesium-rich crystalline silicates absorb stellar radiation poorly ov er visible and 



near-i nfrared wavelengths, where the Sun emits a large fraction of its radiation. IHarker et al. 



(120021 ) argue that Mg-rich pyroxene should have optical con stants very similar to those of 
amorphous Mg-rich pyroxene given by iDorschner et al.l (119951 ) . The imaginary part of the in- 
dex of refraction is as low as ~ 0.0003 at 0.5 fim wavelength for amorphous Mgo.gsFeo.osSiOs, 



so we follow IHarker et al.l (120021) and take this y alue of k for Mg-rich crystalline pyroxene 



which is found in comets (see IHarker et al.ll2002l ). Optical properties for the various silica 
polymorphs at visible and near-infrared wavelengths are not readily available in the liter- 
ature (macroscopic rocks of quartz are often transparent or translucent). We estimate an 
upper limit o n k for both ordinary and extraordinary rays of a-quar tz of 10~^ ( see relevant 
discussion by |Paliklll985l ) and for amorphous Si02 of 10~^ (again see |Paliklll985l ). This sug- 
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gests that, at wavelengths over which the Sun emits most of its power, any sihca in comets 
is either as poorly absorptive as Mg-rich pyroxenes, which are seen in comets, or is even less 
absorptive. Thus, silica may be present in comets but may be at temperatures too low to 
emit enough to be seen in mid-infrared spectra of comet comae. 

A better approach to searching for silica from comets is to analyze physical samples of 
comet dust. In situ measurements of comet dust to date do not determine mineralogy, only 
the abundances of elements constituting comet dust. Such measurements determined the 
rocky material in the coma of Comet Halley to have solar, or chondritic, abundances of the 
major rock-forming elements Mg, Si, Ca, and Fe (jHannerll2003l ). The STARDUST mission to 
Comet 81P/ Wild 2 returned one se veral-micron sized gra in ("Ada") composed of tridymite 
and fayalite (jZolensky et al.ll2006l ). iMikouchi et al.l (120071 ) report a grain of crystalline silica 
that is either tridymite or cristobalite or both from Comet 81P/Wild 2. It remains to be 
seen how much silica will be inferred for this comet based upon analysis of the returned 
samples of its dust. 



As iBradleyl (120031 ) su mmarizes, some Interplaneta ry Dust Particles (IDPs) are believed 
to originate from comets. iRietmeijer fc McKayl (119861 ) reported a 1.5 /im size grain of Si02 
in the phenocryst of the IDP W7010*A2 constituting less than 5% of the abundance of a ll 
minerals in the IDP. This S iQ2 grain is of unknown polymorph (IRietmeijer &: McKaylll985l ). 
Rietmeijer &: McKayl (Il986l ) also reported Si02 fragments of < 10 % abundance in the same 
IDP. These cryst als are subhed ral and ~ 0.1 /im in diameter (IMacKinnon &: Rietmeijer 
19871) . However, iBradlevI (120031 ) does not mention silica as a significant mineral phase in 
IDPs. iGriin et al.l (120011 ). in the discussion of micrometeorite mineralogy and petrography, 
do not list silica as a major phase of micrometeorites, which are IDPs that have landed on 
Earth near the poles. 

In a study o f a lar ge number of Kuiper Belt Objects (KBOs) using near-infrared spectra, 
Barkume et al.l (12008! ) found the spectra were well-fit by models assuming water ice and 
featureless spectral continuum from unknown material. However, other compositions are 
inferred from near-infrared spectra of K BOs, such as metha ne and ethane ice, pyroxenes and 
olivines, and possible organic material (lEmery et al.ll2007l ). It is unknown whether KBOs 
contain silica. 

Cristobalite and tridymite commonly occur in terrestrial siliceous volcanic rocks, such as 
in rapidly-cooled obsidian and rhyolite lava flows of San Juan dis trict, Colorado, USA, b oth 
as the lining of cavities and in the fine-grained groundmass (e.g., iKlein fc Hurlbutll 19771 ). 
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5.3. Origin of Silica in Protoplanetary Disks 



Less than 5% of the Si in the diffuse ISM is in crystalline silicates (ILi &: Drainell200ll ). 
and lKemper et al.l (120041 ) did not detect silica in the interstellar medium 8-12.7 fim absorp- 
tion profile observed towards the Galactic Center. The profile has no 9 or 12.3-12.6 /xm 
features to indicate the presence of amorphous or crystalline silica, and they very nicely fit 



the op tical depth profile without silica. This optical depth profile was shown by lSargent et al. 



(120061 ) to be very similar to the 8-12.7 fim em issivities derived for C oKu Tau/4, DM Tau, and 
GM Aur, whose IRS spectra were shown by lD'Alessio et al.l (120051 ) and lCalvet et al.l (120051 ) 
to be consistent with protoplanetary disks with many-AU-sized almost empty inner holes 
(CoKu Tau/4 and DM Tau) or optically thin inner regions and an almost empty gap (GM 
Aur). The 10 /^m silicate emission from these systems therefore originates from their outer 
disk regions, where the dust is cooler and where dust processing is less likely to occur. Their 
10 fim silicate emission profiles were very well fit with amorphous olivine and amorphous 
pyroxene and negligible silica, similar to the ISM absorption. We infer the average starting 
mixture of grains in protoplanetary disks is like those of CoKu Tau/4, DM Tau, and GM 
Aur and has negligible amounts of silica. Therefore, the silica we see must have formed in 
the disk. 



In the creation of smoke from a laser-ablated enstatite grain, iFabian et al.l (120001 ) found 
that the smallest smoke grains, 10-50 nm in size, were close to Si02 in composition. The 
smallest particles in a smoke created from a forsterite grain were not pure Si02, but had 
Mg/Si ratios between 0.5 and 1. The largest grains in both smokes were close to the com- 
position of their parent enstatite and forsterite grains . This could be an equilib rium process 
related to one in which enstatite melts incongruently (iBowen fc Andersonlll914j ) to forsterite 
plus a high temperature polymorph of silica, either cristobalite or tridymite: 



2MgSi03 ^ Mg2Si04 + Si02 
(Enstatite) (Forsterite) -|- (Cristobalite or Tridymite), 



(3) 



except that in this case amorphous pyroxene, [Mgi_.T;, Fe.r]SiO,s, replaces MgSiO .s, and amor- 
phous olivine, [Mgi_a;,Fea;]2Si04, replaces Mg2Si04 (IBowen &: Schairerl Il935l ). Here, x is 
probably around 0.5 - we will call this reaction "incongruent melting of amorphous py- 
roxene". We invoke amorphous pyroxene as the starting material in this reaction because 
amorphous pyroxene along with amorphous olivine are hypothes ized to be the main ingre- 
dients of the interstellar medium (ISM; see iKemper et al.l |2004| ) and the starting mixture 



of dust for protoplanetary disks (see discussion of unprocessed dust of transitional disks 
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CoKu Tau/4, DM Tau, and GM Aur by ISargent et al.ll2006l ). We also allow that incon- 
gruent melting of enstatite could be happening, where the parent enstatite could itself 
be the daughter product of prior processing of "grandparent" amorphous silicate. Also, 
shocked enstatite grai ns in meteorites are known to have formed cristo balite inclusions 
teenzerara et all l2002h . Shock heating within the protopl anetary disk (11 



Marker fc Desch 



2OO2I ): heating during disk-processing, flares, and lightning (jPilipp et al.lll998l ): heating by 

by th e T 
1996h all 



intense light during the grains' removal from the vicinity of the cent ral starfs) 



Tauri phase solar wind (Ide Pater fc Lissaued I2OOII) or by an X-w i nd (IShu et al. 



may c a use melting of enstatite. E xperiments bv jRietmeiier et al. iHallenbeck et al.l 

fll998l ). Eallenbeck fc NuthI Jl998h . [Rietmeiier et all fll999h . and lRietmeiier et all J2002h also 
found silica grains (amorphous Si02 and tridymite) condensed from vapors of silicate com- 
positions. This suggests silica can form in protoplanetary disks as a smoke or condensate 
from parent silicate dust grains, preferentially of low Mg/Si ratio like pyroxene. However, 
silica is no t expected to condense from vapor of solar composition in chemical equilibrium 
(jGai||200J) in inner disk regions, suggesting the silica grains in the protoplanetary disks in 
our study did not arise from chemical equilibrium condensation in the solar nebula close to 
the Sun. 



5.4. Silica as a Tracer of Dust processing 



Rietmeiier et al. 


(1986 


), 


Hallenbeck et al. 


(1998 




iallenbeck & Nuth 


(1998). 


Rietmeiier et al. 


(1999 


), 


Fabian et al. ( 


2000 


), and 


Rietmeijer et al. 


(2002!) subjected the condensates to hieh 



temperatures to anneal them. The two experiments of annealing at the highest temper- 
atures for the lo ngest durations, t h e hea ting of silica from iron-silicate smokes at 1300K 
for 120 hours by IHallenbeck et al.l (119981 ) and of pure sili ca grains at 1 220K for 5 hours 
bv IPabian eralTj200oh . gave very similar opacity profiles. iPabian et al.l (I2OOOI ) determined 
that their annealed silica was mostly cristobalite with s ome tridymite. The sim ilarity of 
the opacity spectra indicates a similar composition for the IHallenbeck et al.l (119981 ) annealed 
silica. The silica grains remaining after annealing at lower temperature for lesser durations 
in the other annealing experiments of iPabian et al.l (120001 ) remained tridymite and amor- 
phous Si02. Tridymite forms at lower temperatures and does not have an appreciable 16 
/im feature. Therefore the 16 /xm features of IRXS J161410.6-230542 and ROXs 42C, which 
were better fit with annealed silica (which has a 16 /im feature) than with other versions 
of silica, suggests the presence of at least some cristobalite in our silica exemplars. Thus, 
we expect the silica grains in these protoplanetary disks were subje cted to temperatures for 
durations of time like the most extreme annealing experiments of IHallenbeck et al.l (119981 ) 
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and iFabian et al.l fl2000h 



One possibility is that tridymite or cristobalite forms close to the central star, where the 
stellar irradiation heats grains to the high temperatures requir ed for formation of tridymite o r 
cristobalite, and then are transported to cooler disk regions (IBockelee-Morvan et al.ll2002l ). 
Our dust models have annealed silica (cristobalite/tridymite admixture) at temperatures 
ranging from 153K to 1401K, corresponding to distances from the star in the disk of a 
few AU down to the dust sublimation radius (sub-AU range). To explain the existence of 
tridymite or cristobalite in the cool outer disk regions, however, the tridymite or cristobalite, 
once formed, would have to be cooled sufficiently quickly so that their crystalline structure 
does not revert to that of either of the lower-temperature polymorphs of silica, (3- and a- 
quartz, and then transported to the outer disk. It seems unlikely that any of the proposed 
transport mechanisms could move these silica grains to a lower temperature region of the 
disk rapidly enough to quench the cristobalite or tridymite structure in place. 

If the tridymite or crist obalite do not form in the inner disk, this would sugges t transient 
heating events, like shocks (IHarker fc Desclj|2002') or lightning (IPilipp et al.lll998l ). occuring 
in situ in the outer disk. The shock model of IPesch fc Connollvl hom provides cooling sim- 



ilar to that required to explain chondrules - rates of 10-1000 K/hour - which may be quite 
sufficient to allow silica grains to retain their high - temp erature polymorph crystalline struc- 
ture. Furthermore, some of the iDesch &: Connolly J2002h shock models reached temperatures 
high enough that chondrule precursors would be evaporated; this may allow formation of 
silica via condensation. The shock model of production of crystalline silicate grains has been 
invoked to explain the presence of c rysta lline silicates in outflows from Asymp t otic G iant 
Branch (AGE) stars by lEdgar et al.l (120081 ). The lightning theory by lPilipp et al.l (119981 ) also 
attempts to explain chondrule formation and therefore meets similar requirements for cooling 
rates. The heating of pyroxene in some Interplanetary Dust Particles (^I DPs) to tempera- 
tures above 1258K and their subsequent rapid cooling of ~ 1000 K per hour (lAlexander et al. 



20071 ) also suggest heating and rapid cooling in protoplanetary disks, as the chondrules do. 
We note again the chondrules bearing silica by up to 40% by volume and speculate there 
may be a connection between the mechanism that formed submicron silica grains in our 
five silica exemplars and the mechanism that gave rise to the silica-bearing chondrules. An 
experiment like the one that produced the annealed silica whose opacity we use showed that 
such annealed silica cooled from 1223K to 533K over a span of 4 hours 20 minutes with 
average cooling rates dropping from ~ 2000 K per hour initially to ~ 50 K per hour at the 
end (H. Mutschke, private communication). Our data indicate silica grains have experienced 
high temperatures required to form cristobalite or tridymite, and cooled quickly enough to 
quench them. 
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In the various proposed silica-production mechanisms, the sihca forms from parent amor- 
phous sihcate grains of pyroxene ([MgjFeJSiOs) or ohvine ([Mg,Fe] 28104) composition hke 
wha t are inferred to exis t in the ISM from extinction profiles of objects in the Galactic Cen- 
ter (IKemper et al.ll2004j ) or from parent crystalline grains of enstatite (MgSiOs) or forsterite 
(Mg2Si04). Other products arising from the production of silica (Si02) must be rich in 
the leftover magnesium, iron, and/or oxygen. This is seen, for example, in the previously- 
discussed incongruent melting of enstatite or amorphous pyroxene, which produces sihca and 
also forsterite; the forsterite has the magnesium and oxygen from the parent enstatite that 
the daughter silica lacks. This would suggest a correlation between silica abundance and 
silicates with higher stoichiometric ratios of magnesium-to-silicon and oxygen-to-silicon like 
forsterite and the olivines; such correlations should be sought in future dust composition 
studies of protoplanetary disks. 



6. Summary and Conclusions 

We have analyzed the dust composition of five T Tauri stars using their Spitzer Space 
Telescope IRS spectra, with special regard to their silica dust. These five spectra show 
very prominent emission features at 9, 12.6, 20, and sometimes 16 fim wavelength, which is 
characteristic of submicron size silica dust. We have constructed spectral models for each of 
the spectra that include blackbodies at two temperatures, to represent inner and outer disk 
emission from the optically thick disk midplane and blackbody grains, and Planck functions 
at those two temperatures multiplied by scaled dust emissivities, to represent inner and outer 
disk emission from dust in the optically thin disk atmosphere with strong infrared resonances. 
The best fit for a given set of dust opacities happens at the pair of temperatures for which 
the global P^r degree of freedom is minimum. The opacities w e use for non- s ilica d ust 
species are either the same or slightly improved over those used by lSargent et al.l (120061 ). 



We conclude the following: 



• The spectra rule out the highest pressure polymorph of silica, stishovite, and we also 
fairly confidently rule out the intermediate pressure polymorph of silica, coesite. 



Amorphous versions of si lica like amorpho us Si02 (IHenning fc Mutschkd 119971 ) and 
silica-rich glass, obsidian (IKoike et al.lll989l ) are quite firmly ruled out to explain the 
silica features in our five silica exemplars due to their inability to fit the 12.6 fim 
features seen in these spectra. The amorphous Si02 features peak at too short of a 
wavelength for this feature, are too wide, and are too weak with respect to the 9 /xm 
feature. 
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• The opacity of anne aled silica obtained by heating amorphous Si02 for 5 hours at 1220K 
(IFabian et al.ll2000l ) prov ides the best fit to the spectra of all five of the strong silica 
exemplars in this study. iFabian et al.l (|2000| ) report this annealed silica to be mostly 
cristobalite but partly tridymite, so these silica polymorphs are the best candidates 
to account for the silica features seen in t he spectra of our five sili ca exemplars. This 
opacity is very similar to that reported by iHallenbeck et al.l (119981 ) for silica annealed 
at 1300K for 120 hours. 

• The opacity of a-quartz at intermediate temperatures (505K, 620K) provided a good 
fit to the 12.6 /im feature of ROXs 42C, but the fits to features at other wavelengths 
characteristic of silica (9, 16, 20 /im) were not as good as the fits obtained from using 
annealed silica, indicating a-quartz is not the dominant silica in these disks. /5-quartz, 
however, can be fairly firmly ruled out, as its feature at 12.7 fim cannot fit the 12.6 
fim feature we observe in our dust excess spectra, regardless of the assumed shape 
distribution of the /9-quartz grains. 

• Neither amorphous nor crystalline silica are seen in interstellar medium dust, so the 
silica present in protoplanetary disks of the five silica exemplars analyzed in this study 
must have formed in the disks. 

• The silica grains giving rise to emission from our five silica exemplars may have formed 
in similar conditions to those that gave rise to Ada, a grain composed mostly of 
tridymite surrounded by fayalite, found in the coma of Comet Wild/2 and the sil- 
ica polymorphs found in primitive meteorites. 

• Tridymite and cristobalite are common products of annealing of amorphous dust 



smoke s and condensa t es in l a boratory experiments (e. g, iRietmeijer et al.l (119861 ) , iRietmeijer et al 
(119991 ) ■ IFabian et al.l (120001 ) . IRietmeijer et al.l (120021 )) designed to mimic protoplane- 
tary disk conditions thought to produce crystalline silicates by annealing amorphous 
silicate smokes. Cristobalite, the major component of the annealed silica we use, is 
present in silica grains annealed at temperatures above 1220K for durations more than 
5 hours and has a 16 /im feature in its opacity profile that is needed to fit the spectra 
of our silica exemplars. Silica, as a tridymite or cristobalite polymorph, might form in 
protoplanetary disks in manners similar to these laboratory experiments. Protoplane- 
tary disk temperature and heating duration conditions giving rise to the silica we see 
in our sample of protoplanetary disks must be more like the most extreme conditions 
from the annealing experiments giving rise to cristobalite. 



• The silica grains may form via the incongruent melting reaction of amorphous pyroxene 
or the incongruent melting reaction of enstatite. 
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• The tridymite or cristobalite grains, once formed, must be cooled quickly enough to 
retain their crystalline structure and not revert to the lower temperature polymorphs 
of silica, (3- and a-quartz. Mechanisms like spiral shocks and lightning devised to 
explain the rapidity with which chondrules cool may also satisfy the rapidity with 
which tridymite or cristobalite must cool in order to retain its crystalline structure. 
This would support the idea of in situ formation of the annealed silica. 
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Appendix 

For ZZ Tau, average flux densities for each DCE of each order of each nod were com- 
puted, averaging over 5.18-6.76 fim wavelength for SL second order, 7.57-8.37 /im wavelength 
for SL bonus order, 7.57-11.81 fim wavelength for SL first order, 14.15-17.97 fim wavelength 
for LL second order, 19.82-21.19 /im wavelength for LL bonus order, and 22.02-29.99 /im 
wavelength for LL first order. Each DCE of ZZ Tau was then scaled so that its average fiux 
was the same average fiux as the other DCEs of the same order. For SL second order, all 
DCEs were scaled so that their average fiux was the same as that of the second nod from 
the Campaign 4 AOR; for SL bonus order, the average fiux to match was that of the first 
nod of the Campaign 4 AOR; for SL first order, also the first nod of the Campaign 4 AOR; 
for LL second order, the second DCE of the second nod of the Campaign 29 AOR; for LL 
bonus order, the first DCE of the first nod of the Campaign 29 AOR; and for LL first order. 
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the second DCE of the first nod of the Campaign 29 AOR. In almost all cases, the average 
fluxes of these named DCEs were the highest of their given order; in the case of LL second 
order, the average flux of the named DCE was only 1% discrepant from the highest average 
flux of all LL second order DCE spectra. The range of scalars applied to each DCE of ZZ 
Tau is provided also in Table 3. 
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Fig. 1. — (Color on-line only) Comparison of op acities of a-quartz in the CDE and CDE2 
shape distributions using optical properties from iGervais fc Pirioij (119751 ) for a-quartz at 
295K. Note the shifting of the long-wavelength sides of the 9 and 20 fim features, but not of 
the weaker 12.6 and 14 /xm features, to longer wavelengths. 
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Fig. 2. — (Color on-line only) Comparison between opacity profiles computed for grains 
in CDE shape distribution using optical constants of a-quartz and /3-quartz from Gervais 
and Piriou 1975. The opacity profile of a-quartz at 295K is plotted along with the average 
of the opacity profiles of a-quartz at 505K and 620K, the average of the opacity profiles of 
a-quartz at 785K and 825K, and the average of the opacity profiles of /3-quartz at 860K and 
975K. 
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Fig. 3. — (Color on-line only) Comparison of opacities of obsidian and amorphous Si02 in 
the CDE shape distribution and annealed silica. Note the ratio of maximum opacity in the 
12.3 fim feature to the maximum opacity in the 9 /im feature for obsidian and amorphous 
Si02 is lower than that of the 12.6 /im feature to the 9 /im feature for annealed silica. 
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Fig. 4. — (Color on-line only) KBr effect on /3-quartz at 975K. The peak wavelengths 
of the 9 /im features of /?-quartz grains in both shape distribution in KBr and of annealed 
silica grains (in KBr) are very close. The CDE2 shape distribution of /5-quartz in vacuum 
gives a 9 /im peak ~ 0.2 fim shortward of the same grains in KBr, while the same peak for 
/9-quartz in the CDE in vacuum is not much different in central wavelength from the peaks 
of /3-quartz in KBr. The ~ 20 fim complexes of the various /5-quartz opacity curves are 
similarly affected; the weaker 12.7 /^m feature is hardly affected in terms of peak width and 
central wavelength by shape or KBr/vacuum considerations at all. 
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Fig. 5. — (Color on-line only) Model fit to ROXs 42C using annealed silica, minimizing 
chi-squared from 6 to 37 fim, comparing fit including forsterite in the tCDE (black solid 
line; see text for description of tCDE) and CDE (orange dashed line) shape distributions. 
The model was constructed as described in Section 3 using parameters listed in Table 2. 
All non-silica model components for the model using forsterite grains in the tCDE shape 
distribution are plotted as solid lines at bottom of graph. Model components from hot dust 
are plotted in blue, and those from cold dust are plotted in red. Silica components are 
plotted with a dashed line for the hot silica and a dotted line for the cold silica. Note the 
lack of the 6.2 fim PAH feature. 
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Fig. 6. — (Color on-line only) Model fits to ROXs 42C. The models are solid black lines, 
each model using a different kind of silica and paired with the spectrum of ROXs 42C. Each 
model/spectrum pair is translated vertically by 0.2 Janskys. The models use, from bottom 
to top in the plot: obsidian, amorphous Si02, quartz in the CDE2 shape distribution, quartz 
in the CDE shape distribution, and annealed silica. 
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Fig. 7. — (Color on-line only) Comparison of chi-squared per degree of freedom for dust 
emission models using different types of silica for the spectra of the five objects modeled, 
including the dereddened spectra obtain by two different computations of optical depth at 
9.7 /im of the material in front of IRS 49. The type of quartz - /3-quartz and an average of 
a-quartz at 785K and 825K - is denoted by /5 and an respectively. 
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Fig. 8.— (Color on-line only) Best fits to ZZ Tau (top) and IRXS J161410.6-230542 
(bottom) using annealed silica. Same color and linestyle scheme for model components as 
for Figure 5. The dominant non-continuum dust component in the 8-12 fim region of ZZ 
Tau is large amorphous pyroxene. 
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Fig. 9. — (Color on-line only) Best fits to IRS 49 (top) and T51 (bottom) using annealed 
silica. Same color and linestyle scheme for model components as for Figure 5, except for 
the orange points being the spectrum of IRS 49 extinction corrected assuming Ay/rg j = 25 
instead of 18 and the solid violet line running through the orange points being the model fit 
using annealed silica to this differently-corrected data. 
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Table 1: Observations and Stellar Characteristics 









R.A. 


Dccl. 


Spectral 












Object 


AOR 


# in AOR 


(J2000) 


{J2000) 


Typo 


Av 




^star 


dist. 


Rcfs. 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 



IRXS J161410. 6-230542 5206272 ■■■ 16 14 11.08 -23 05 36.2 KO 1.48 4963 5.18 144 1,2,3 

IRS 49 12698624 6 16 27 38.31 -24 36 58.7 K5.5 10.7 ■ ■ ■ ■ . ■ 140 4,5,6 

ROXs 420 12676352 1 16 31 15.75 -24 34 02.2 K6 1.9 4205 8.54 140 6,7 

T51 12696576 4 11 12 24.41 -76 37 06.4 K3.5 ■ ■ ■ ■ ■ ■ 160 8,9 

ZZ Tau 3546880,16261376 ■_ 04 30 51.38 +24 42 22.3 M3 1.44 3470 4.90 140 10,11,12 



Note. — Col. (1): Object name. Col. (2): Spiteer Astronomical Observation Request number. Col. (3): 
If multiple objects observed in AOR, then specifies which one observed. If only one object in AOR, then no 
number specified. Col. (4): Right Ascension of object in J2000 coordinates. Col (5): Declination of object 
in J2000 coordinates. Col. (6): Spectral type of object. Col. (7): Extinction at V band in magnitudes. 
Col. (8): If stellar photosphere emission was subtracted from the IRS spectrum, then this entry is the 
effective temperature of the star which is the temperature of the blackbody representing stellar photosphere 
emission normalized to the dereddened J-band flux density. If stellar photosphere emission not subtracted, 
then no effective temperature specified. Col. (9): If stellar photosphere emission subtracted, then solid 
angle in units of 10~^^ steradians of blackbody representing stellar photosphere emission. Obtained by 
dividing dereddened J band flux density by Planck function evaluated at isophotal wavelength of J band 
and at temperature specified in previous column; if no stellar photosphere emission subtracted, then no 
solid angle specified. Col. (10): Assumed distance t o object, in par s ecs. C ol. ( 11): References from w hich 



spectral type, Ay, an d dis tance are obta i ned. 1 - iPreibisch et al.l (119981) : 2 - iPascucci et al.l (120071) : 3 



deZeeuwetal 



( 19991): 4 - IWilking et al.l (l2005l): 5 - iMcClure eTall (l2008l) : 6 - iBontemns et all toOlh : 



Bouvier fc Appenzelleil (Il992l): 8 - iLuhmanI j20oi): 9 - iLuhmanI (I2OO8I): 10 - iHerbig fc Belli (Il988l): 11 



7 



Furlan et al 



(|2006l ): 12 - iKenvon et all |l994l ). For I RXS J161410.6-230542 Tef f comes fromlPascioccLelal 
(|2007[ ): for ROXs 42C and ZZ Tau, Teff comes from IKenvon fc HartmannI (jl995h based on the spectral type 
from column 4. 



Table 2. Dust Model Parameters 













small 


small 


large 


large 






Total 






Temp 




silica 




Amorphous 


Amorph. 


Amorph. 


Amorph. 


Orystalline 


Crystalline 


Dust 






\'^) 


Odd 


type 






Ill T\^i n 




V-/11 V lllC 


JT V 1 LIACllC 


a IfcJl 1 It) 


iVlcloo 


d.o.f. 


{'■) 








\°) 




\') 


C8'l 










IRXS J161410.6 


305 


57.5 ± 0.8 


annSil 


35.7 ± 17.3 





12.0 ± 16.3 








26.1 ± 14.5 


26.3 ± 14.2 


0.815 


1.1 


-230542 


1401 


0.329 ± 0.020 




5.6 ± 3.2 








94.4 ± 16.0 











0.0474 




IRS 49'' 


153 


540 ± 34 


annSil 


6.2 ± 2.5 





7.8 ± 3.6 





70.6 ± 6.4 


12.2 ± 2.6 


3.2 ± 2.5 


176 


5.5 




751 


3.36 ± 0.11 




14.0 ± 1.5 








86.0 ± 4.5 











0.904 




IRS 49' 


160 


404 ± 26 


annSil 


6.6 ± 2.3 











79.7 ± 6.2 


10.1 ± 2.3 


3.6 ± 2.2 


150 


6.3 




796 


3.02 ± 0.09 




16.2 ± 1.9 








83.8 ± 5.2 











0.569 ' 
25.3 




ROXs 42C 


191 


429 ± 9 


annSil 


20.0 ± 5.4 


11.6 ± 7.7 


32.1 ± 7.4 








14.6 ± 5.1 


21.8 ± 5.5 


3.2 




1051 


0.662 ± 0.046 




13.8 ± 2.6 


5.4 ± 3.8 


16.5 ± 4.5 


57.1 ± 6.8 





5.8 ± 3.3 


1.5 ± 3.2 


0.172 1 




ROXs 42C 


172 


613 ± 16 


oqCDElt 


9.5 ± 6.4 





37.7 ± 10.1 








42.9 ± 10.5 


9.9 ± 5.9 


32.9 


5.0 




601 


4.18 ± 0.15 


oqCDEht 


18.4 ± 3.8 





69.1 ± 11.4 


0.1 ± 9.5 








12.3 ± 4.4 


0.393 




ROXs 420 


191 


421 ± 9 


aqCDE21t 


15.1 ± 4.4 





15.1 ± 5.5 





39.5 ± 6.8 


16.3 ± 4.2 


14.1 ± 4.0 


30.3 


4.5 




1151 


0.707 ± 0.039 


/3qODE2 


3.9 ± 2.3 


70.3 ± 13.9 











19.4 it 5.9 


6.4 ± 5.1 


0.0900 




ROXs 420 


191 


400 ± 9 


amsil 


17.9 ± 5.1 





33.8 ± 7.0 








29.7 ± 6.0 


18.6 ± 4.8 


27.4 


5.2 




751 


1.86 ± 0.09 




10.3 ± 3.3 


59.0 ± 12.4 


17.4 ± 7.8 








8.2 ± 5.8 


5.2 ± 5.3 


0.193 




ROXs 420 


210 


275 ± 6 


obsid 


25.7 ± 4.1 


21.8 ± 4.9 


22.5 ± 4.5 








15.7 ± 3.4 


14.3 ± 3.2 


24.8 


5.4 




1401 


0.376 ± 0.028 




25.1 ± 8.1 





62.8 ± 18.1 








9.9 ± 8.1 


2.3 ± 7.8 


0.0427 




T51 


267 


39.2 ± 1.7 


annSil 


4.8 ± 3.6 











70.9 ± 9.3 


15.2 ± 3.9 


9.2 ± 3.6 


6.14 


3.6 




1351 


0.792 ± 0.037 




9.1 ± 1.1 








88.2 ± 4.0 





1.8 ± 1.5 


1.0 ± 1.4 


0.356 





Table 2 — Continued 



small small large large Total 

Temp silica Amorphous Amorph. Amorph. Amorph. Crystalline Crystalline Dust 

Object (K) f^ss type Silica^ Pyroxene'' Olivine'^ Pyroxene'^ Olivine^ Pyroxene^ forsterite^ Mass gT^T' 

(1) (2) {3^ (4) (5) (7^ {9^ 00^ (11) (12) (13)' 

ZZ Tau 210 62.9 ± 1.5 annSil 45.8 ± 23.8 54.2 ± 27.1 1.46 3.3 

1251 0.123 ± 0.007 "" 4.2 ± 1.9 5.5 ± 3.2 85.5 ± 8.2 4.2 ± 2.9 0.5 ± 2.8 0.0290 



Note. — Col. (1): Object name. Col. (2): One of two dust model temperatures (Kelvin) . Col. (3): Solid angle, fisB, of blackbody of temperature specified 
in Col. (2) representing continuum emission, expressed in units of 10~^^ steradians. Col. (4): Type of silica used in model. Cols. (5)-(ll): mass percentages 
of all dust mass at temperature specified in Col. (2). One dust model is completely specified by two adjacent rows - the row following the object's name 
and the row beneath that one. Col. (12): Total dust mass at one temperature in 10"'* lunar masses, computed assuming distances to each object as listed 
in Table 1. Col. (13): per degree of freedom, determined over 7.7 < A < 23 /im. 

'^Optical properties and opacities for various types of silica from references given i n text. In the table, "an nSil" means annealed silica, "aqCDIflt" means 
an average of the CDE opacity curves obtained for a-quartz at 505K and 620K from lCervais fc Piriou ( 1975 ): " aqCp E21t" means the same, but^r CDE2; 
"aqCDEht" means an average of CDE opacity curves for a-quartz at 785K and 825K from Gervais fc Piriou ( 1975h : "/3qCDE2" means the saqie, but for 
CDE2; "obsid" means obsidian; and "amsil" means amorphous Si02. 

''Optical constants for amorphous pyroxene Mgo.yFco.aSiOs from Dorschner et al. ( 19951 ). assuming CDE2 ( Fabian et al. 2001 ) 

Optical constants for amorphous olivine MgFeSi04 from IPorschner et al. ( 19951 ). assuming CDE2 

''Optical constants for amorphous pyroxene Mgo.yFco.aSiOs from lDorschner et al.l (|l995| ). using the Bruggeman EMT and Mie theory (Bohren fc Huffman 
1983) with a volume fraction of vacuum of f — 0.6 for porous spherical grains of radius 5 /im 

■^Optical constants for amorphous olivine MgFeSi04 froml Dorschner et al. ( 1995 ). using the Bruggeman EMT and Mie theory (Bohren & Huffman 1983) 
with a volume fraction of vacuum of f — 0.6 for porous spherical grains of radius 5 /zm 

'Opacities for crystalline pyroxene Mgo.gFeo.iSiOs from lchihara et aL ( 2002 ) 

^Optical constants for 3 crystallographic axes of forsterite, Mg2Si04, from Sogawa et al. ( 2006h . assuming tCDE (see discussion in Section 3.3 of this 
paper) shape distribution. 

''Extinction corrected using Ay /rg. 7=18; r9.7=0.6 



'Extinction corrected using Ay /rg. 7=25; tq,7=0A 



I 

I 
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Table 3: Mispointing Corrections 





Correction to 
















Object 


SL2nodl 


SL2nod2 


SLlnodl 


SLlnod2 


LL2nodl 


LL2nod2 


LLlnodl 


LLlnod2 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


IRXS J161410.6-230B42 


1.03 


1.03 


1.03 


1.03 


1.00 


1.00 


1.00 


1.00 


IRS 49 


1.00 


1.04 


1.00 


1.04 


1.00 


1.04 


1.00 


1.01 


ROXs 42C 


1.00 


1.00 


1.00 


1.00 


1.00 


1.02 


1.00 


1.00 


TBI 


1.00 


1.00 


1.00 


1.00 


1.00 


1.01 


1.00 


1.01 


ZZ Tau 


1.00-1.10 


1.00-1.10 


1.00-1.14 


1.02-1.07 


0.99-1.09 


1.00-1.09 


1.00-1.08 


1.02-1.08 



Note. — Col. (1): Object name. Cols. (2)-(9): Multiplicative scalars applied to the spectrum of one order 
of one nod for a given object to match the flux density of the other nod as described in the text. SL2 means 
Short-Low order 2, SLl is Short-Low order 1, LL2 is Long-Low order 2, LLl is Long-Low order 1. For ZZ 
Tau, the range of scalars used for each order of each nod, derived as described in the text, is provided. 



